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Potentials between static colour sources in eight differ- 
ent representations are computed in four dimensional SU{3) 
gauge theory. The simulations have been performed with the 
Wilson action on anisotropic lattices where the renormalised 
anisotropics have been determined non-perturbatively. After 
an extrapolation to the continuum limit we are able to ex- 
clude any violations of the Casimir scaling hypothesis that 
exceed 5 % for source separations of up to 1 fm. 

PACS numbers: 11.15.Ha, 12.38.Gc, 12.38.Aw, 12.39.Pn 



I. INTRODUCTION 

Non-perturbative QCD effects in general and the na- 
ture of the confinement mechanism in particular are the- 
oretically challenging. At the same time these aspects 
are important for high energy and low energy particle 
and nuclear phenomenology. Several models of non- 
perturbative QCD have been proposed whose predictions 
happen to differ from each other substantially in some 
cases. Prominent examples are bag models strong 
coupling and flux tube models |^-^, bosonic string mod- 
els pO^O|, the stochastic vacuum model [121-^^, dual 
QCDl|l|^|l3, the AbeUan Higgs model ||l|]l9l|7and in- 
stanton based models p0|-|2^ . Lattice simulations of in- 
teractions between static colour sources offer an ideal en- 
vironment for discriminating between different models of 
low energy QCD and to learn more about the confine- 
ment mechanism. They are easily accessible analytically 
and at the same time very accurate Monte Carlo predic- 
tions can be obtained |p3| , p^ . 

Despite the availability of a wealth of information on 
fundamental potentials, only few lattice investigations 
of forces between sources in higher representations of 
SU{N) gauge groups exist. Most of these studies have 



been performed in SU{2) gauge theory in three 25 - 2q] 
and four |29|-|3^ space-time dimensions. Zero temper- 
ature results for four dimensional SU (3) can be found 
in Refs. |37-4^] while determinations of Polyakov line 
correlators in non-fundamental representation have been 
performed at finite temperature by Bernard ps|j30[ | for 
SU{2) and in Refs. for Suli) gauge theory. 



In our study we shall see that the so-called Casimir 
scaling hypothesis is rather accurately represented 
by the lattice data while models predicting a different be- 
haviour are definitely ruled out. Casimir scaling means 
that potentials between sources in different representa- 
tions are proportional to each other with their ratios 
given by the respective ratios of the eigenvalues of the 
corresponding quadratic Casimir operators, which is ex- 
act in the case of two dimensional Yang-Mills theories. 
Our result is of particular interest with respect to recent 
discussions of the confinement scenario [^^-^ . 

At distances r > rc « 1.2 fm non- fundamental 
sources will be screened and "string breaking" effects will 
be encountered that are incompatible with Casimir scal- 
ing. In the present study we restrict ourselves to dis- 
tances smaller than the string breaking scale Vc- 

This article is organised as follows: in Sec. || the lattice 
methods that we apply and our notations are introduced. 
A determination of the renormalised anisotropies and lat- 
tice spacings is presented in Sec. HI. The potentials are 
then determined in Sec. |^ before we conclude with a 
brief discussion. 



II. NOTATIONS AND METHODS 

We denote the energy of colour sources, separated by a 
distance r, in a representation D = 3, 6, 8, 10, ... of the 
5C/(3) gauge group by Voir,^), where /i denotes some 
cut-off scale on the gluon momenta, for instance an in- 
verse lattice spacing, /i = it /a. We shall also use the 
subscript "F" to label the fundamental (3) representa- 
tion (or we may just omit the subscript in this case). 

TABLE I. Group theoretical factors for SU{'i). D is the 
dimension of the representation, (p, q) are the weight factors, 
z — exp(27ri/3), and dn ~ Cd/C'f denote ratios of quadratic 
Casimir charges. 
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The static potential, 

VD{r,^Ji) Voir) + Fd.sciKm), 



(1) 



can be factorised into an interaction part, VD(r), and a 
self energy contribution, VD,seif(M) that will diverge like 
/i/ In /i as /i — *■ cxD while Vd('') will assume universal val- 
ues. 

A (dimensionless) lattice potential Vd (R-, a) will re- 
semble the corresponding continuum potential up to lat- 
tice artefacts, 



VoiRa) = fl-i \ VD{B.,a) - Vb,scif(a) 
1 + fD{Ra,fC)a'' 



(2) 



where is a positive integer number that will in general 
depend on the lattice action employed. We are concerned 
with Wilson-type gluonic actions |Q. In this case, v = 2. 
Note that the coefficient function / only depends on the 
combination r = Ra and on the direction of R but not 
on R itself. This guarantees that lattice artefacts are 
reduced as r ^ a is increased. 

The static potentials are obtained from fits to 
smeared [^,^^ Wilson loops for T > T^^in where T^in 
depends on R, the statistical errors of the Wilson loops 
and the smearing algorithm employed. We define a Wil- 
son loop in representation D, 



Wd(R,T) =Tr 



n ^D,n.f^ I ' (3) 

'.,n)e5CCR,T) 



where (5C(R, T) denotes the oriented boundary of a (gen- 
eralised) rectangle with spatial extent R and a temporal 
separation of T lattice spacings. (n, /i) denotes an ori- 
ented link connecting the site n with n + fi, n is an in- 
teger four-vector that labels a lattice site and /i is a unit 
vector pointing into a direction, G {1,2,3,4}. "Tr" is 



the normalised trace, TyuIu 



-tr l£) = 1, Nd is the 



dimension of the representation D and, 



Ur 



V |exp 



a(n+fi} 



(4) 



denotes a link variable in representation D, where 
is a generator of the gauge group. Our conventions are 
[T^,T^] = ifabcT^, where fabc are totally antisymmet- 
ric real structure constants. The normalisation is such 
that tr T^T^ = Sab/^. Now: 

{Wd{R, T)) = cz5(R, a) exp [-Vi)(R, a)T^ (T -> oo). 

(5) 

The use of smeared Wilson loops turns out to be more 
suitable for numerical simulations than implementing the 
definition of Eq. (||) ; the spatial pieces of the Wilson loop 



are replaced by linear combinations of various paths that 
models the ground state wave function. As a result the 
overlap of the creation operator with this ground state, 
CD, is enhanced and the T — *■ cx3 limit can effectively be 
realised at moderate T values. 

In numerical simulations one observes that the statis- 
tical error AW{T) of the expectation value of a smeared 
Wilson loop {W{T)) only weakly varies with T [|^. From 
Eqs. (|l|) and (||) we therefore obtain the relation, 

AW{T) 



{Win 



oc exp[y(a)T] cx exp Vsou{a)T (T oo) 



(6) 



for the relative errors (which are directly proportional 
to the statistical uncertainty of the potential values). In 
tree level perturbation theory one finds. 



VD,scu{a) = cCDgHa), c = 0.252731 



(7) 



the self energy is proportional to the eigenvalue of the 
quadratic Casimir operator Cd = Ti dT^T^ of the rep- 
resentation. This means that statistical errors will in- 
crease significantly as we investigate higher representa- 
tions of the sources with bigger Casimir charges. 

This self energy related problem motivates us to intro- 
duce an anisotropy parameter ^ = aa/dr ~ 4 between 
spatial lattice resolution Ua- and temporal lattice con- 
stant Gr- This results in a reduction of the self energy, 
Vd.scU = cCd £,~^g^, and therefore of the relative errors 
of smeared Wilson loops. However, at the end of the day 
we wish to measure distance and potential in the same 
units flo-. This means that we do not gain anything from 
the factor ^ within the above expression. Nevertheless, a 
^ > 1 still results in a reduced effective at fixed Ua- ■ Of 
course, equally well we could just have increased the size 
of our statistical ensemble by a factor four and worked 
at ^ = 1. 

Our main motivation for introducing an anisotropy is 
the possibility of reducing lattice artefacts. These are 
most prominent at small and at large distances: as long 
as r is not much larger than Oo- the cubic lattice struc- 
ture is clearly visible. In lowest order perturbation theory 
these violations of rotational symmetry only depend on 
R and while the order 5^ coefficients exhibit a weak 
dependence on ^ too ||5^. While we cannot hope to sig- 
nificantly reduce these small distance effects without de- 
creasing flcr , it is clear that on a lattice with temporal res- 
olution Qr one cannot reliably resolve masses m ^ o-t^- 
However, at r 3> and in particular for representations 
D with large Casimir charges situations, VD{r)aa 3> 1, 
are easily encountered, unless ^ S> 1. 

Introducing an anisotropy also means that within any 
physical t window we have more data points at our dis- 
posal. While this might help to gain more confidence in 
identifying effective mass plateaus we find that the ad- 
ditional data points are highly correlated and add little 
extra information, at least when one is only interested in 
the mass of the ground state within a given channel. 



2 



Our action reads, 



S 



^^Tr[/„,„ +eo^TrC/„ 



(8) 



where /? — 2N/[g'^{acr, £,)] is defined through the lat- 
tice couphng and j,j e {1,2,3}. Un,fj,i^ — 
Un,fj.Un+ii,iyUl^^- dcuotes the product of four hnk 
variables around an elementary square, the plaquette. 
With the above anisotropic Wilson action, the lead- 
ing order lattice artefacts are proportional to and 
arCia = OctC"^- This means that along a trajectory 
of constant ^ the continuum limit will be approached 
quadratically in Oo-. The relationship between the bare 
anisotropy appearing in the action Eq. (^) and the 
renormalised anisotropy ^ is known in one loop pertur- 
bation theory ^ = ^o[l + ci(^o)5^ H ]■ In Sec. |ll 

we will discuss our non-perturbative evaluation of ^. Of 
course the function $(17^,^0) is not unique and different 
non-perturbative definitions might differ from each other 
by terms of order a„ . 

Perturbation theory yields the relation between poten- 
tials in different representations D, 



(9) 



where c?d = Cd/Cf- Table || contains all representa- 
tions D, the corresponding weights {p, q) and the ratios 
of Casimir factors, do, for p + q < 4:. In SU{3) we have 
Cf — 4/3, and z = exp(27ri/3) denotes a third root of 
1. Eq. (||) is known to hold to (at least) one loop (order 
g^) perturbation theory at finite lattice spacing p7[| and 
two loops (order g^) in the continuum limit of four 
dimensional Yang-Mills theories. The main purpose of 
this article is to investigate non-perturbatively to what 
extent the "Casimir scaling" relation Eq. (||) is violated. 

In what follows U denotes a group element in the fun- 
damental representation of SU{3), for instance the prod- 
uct of link variables around a closed contour. The traces 
of Ud in various representations, Vd — tr Ud, can easily 
be expressed in terms of traces of powers of U, 



V3 



trU, 

^3 



(I^3P 



1 



- [{trUf + trU'] 



V6 

Vi5a=tTU*Ve~tTU., 

Vio 

V24 = trU*Vio-Ve, 
V27 



I [(tr Uf + 3tiUtiU^ + 2 tr U^] 
6 



— [(tr + 6(tr Uftr + 3(tr 
8trC/trC/^ + 6trC/^] . 



(10) 
(11) 
(12) 
(13) 
(14) 

(15) 
(16) 

(17) 



Note that TidUd = -jq^trUD = Vd/Nd- Hence, the 
normalisation of Vd differs by a factor No from that of 



the Wilson loop Wd of Eq. (_ 
U ^ zU, Wd transforms like, Wd z- 
sentations with z^^"^ = 1 have zero triality. 



Under the replacement, 
p-'^Wd- Repre- 



III. DETERMINING ANISOTROPY AND 
LATTICE SPACING 

We simulate SU{3) gauge theory at the parameter val- 
ues (/3,Co) = (5.8,3.10), (6.0,3.20) and (6.2,3.25). From 
exploratory simulations with limited statistics and the 
published data of Refs. p^-|61| we expect to find renor- 
malised anisotropics ^ « 4 at these combinations. At the 
above parameter values volumes of L^. x Lr — 8'^ x 48, 
12'^ X 72 and 16'^ x 96 lattice sites have been realised, 
respectively. These volumes were chosen to keep the lat- 
tice extent about constant in physical units. In addition 
a volume of 12"^ x 48 lattice sites has been simulated at 
(3 = 5.8 to investigate possible finite size effects. 

The gauge configurations have been obtained by ran- 
domly mixing Cabibbo-Marinari style |^ Fabricius- 
Haan heatbath [^ and Creutz overrelaxation ]6^ ] 
sweeps, where we cycled over the three diagonal SU{2) 
subgroups. The probability of a heatbath sweep was set 
to be 1/5. During each sweep the sites were visited sub- 
sequently for each of the four space-time directions of 
the links in sequential order. Measurements were taken 
after 2000 initial heatbath sweeps and the gauge configu- 
rations are separated from each other by 200 sweeps. In 
doing so, we did not find any signs of autocorrelation or 
thermalization effects for any of the investigated observ- 
ables at any of the simulated parameter sets. In the case 
of /3 = 6.0 one set of configurations was generated on a 
Sparc station after a cold start while another set of con- 
figurations was generated on a Cray J90, starting from 
a hot, random configuration. No statistically significant 
deviations between these two data sets were found either. 
We display our simulation parameters in Table |l[ riconf 
denotes the number of statistically independent config- 
urations analysed in each case while rp ~ 0.5 fm is the 
Sommer scale parameter ||65|, implicitly defined through 
the static potential. 



dV{r) 



dr 



= 1.65. 



(18) 



TABLE 11. Simulation parameters, lattice spacings and lin- 
ear lattice extents. 



/3 


L„ 


^COIlf 


ro/aa- 


ro/ar 




5.8 


8 


633 


3.074(29) 


12.45(11) 


2.60(3) 


5.8 


12 


139 


3.052(34) 


12.52(13) 


3.93(4) 


6.0 


12 


159 


4.437(71) 


18.12(24) 


2.70(5) 


6.2 


16 


133 


6.106(66) 


24.16(23) 


2.62(3) 
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FIG. 1. Ground state overlaps for the case of the funda- 
mental potential determined from Wilson loops in three dif- 
ferent orientations at /3 = 6.2, = 3.25. 

We label quantities associated with the fine grained 
direction by an index r while a refers to coarse grained 
directions. On an anisotropic lattice various ways of as- 
sociating the sides of smeared Wilson loops with these 
directions exist: Wra^r / a-nt/ aa), Waa{'>^ I ^a^t/ a^) and 
W„r{Y/acr,t/ar). In the case of Wra as well as for Waa 
the time coordinate points into a a direction. The spa- 
tial coordinate is identified with the r direction in the 
first case and a a direction in the latter case. While the 
spatial connections within Wra are parallel to the r axis, 
in the case of Wo-o- we realise planar off-axis configura- 
tions R II (L 1, 0) and R || (2, 1, 0), in addition to on-axis 
separations^]. Finally, within W^r the time coordinate is 
taken along the fine grained dimension and the spatial 
coordinate coarse grained. We determine W^t for the 
standard separations |||, R || (1,0,0), (1,1,0), (2,1,0), 
(1,1,1), (2,1,1), (2,2,1). 

In the spatially isotropic situation (War) we iteratively 
construct fat links in the standard way by replacing a 
given link by the sum of itself and the neighbouring four 
spatial staples with some weight parameter, a > 1, 



Psuis) ct 



(19) 



Un,jUn+j,iU^^-j + U\_~JJn-'i,iUn-'i+j,i 



Psu(3} denotes a projection operator, back onto the 
SU{3) manifold. We employ the definition j2^, U = 
PsuisA-A) G SU{3), ReTrL/At = max and iterate 
Eq. ()l|) 26 times with a = 2.3. 



We use a somewhat different novel smearing algorithm 
in the case of Wra and Waa where the spatial volume 
is anisotropic with one fine and two coarse directions: 
when one only considers links parallel to the one being 
replaced, Eq. (jl^) resembles a two dimensional diffusion 
process: U Psu{3) [(q;+4+V2)C/]. We are interested to 
maintain an isotropic propagation of the link fields when 
an anisotropy parameter ^ > 1 is introduced. Following 
the above diffusion model this is achieved by replacing 
Eq. (|l|) with. 



PSU{3) (a Un^^ + Fn^] + S.^Fn,r) 



Un,T * P; 



SU(3) 



ia + 2e-2)Un,r + J2^"'^ 



(20) 
(21) 



where i,j e 7^ We perform 22 iterations of 

Eqs. ( p0| ) - ( pi| ) with a = 3.7. Indeed, in doing so we 
find very similar overlaps with the physical ground state, 
C(j(j(R) « Cro-(R), where c = c^? S [0, 1] is defined in 
Eq. (|) and WraCR^T - 0) = W,.(R,r = 0) = 1. 
This is illustrated in the comparison of data obtained at 
P = 6.2, ^0 = 3.25 of Fig. |l|. We have not been able, 
however, to sustain the high overlaps achieved for War 
(triangles) for Wra or Waa- The situation at the other 
(/3,^o) combinations is similar. 

From the asymptotic behaviour of the different Wil- 
son loops at large temporal separations t, three lattice 
potentials can be determined: 



a'^'^Vrair/ar) 

a^^Vaa{Y/aa) 
a~'^Var{r I aa) 



Mm —\nWra{r/ar,t/aa), (22) 

t^oo at 

Mia ^\nWaa{rlaa,t/aa), (23) 

t^oo dt 

Mm ^\nWar{r/aa,t/ar). (24) 



Note that while Vra and Vaa are given in units of Ua , Var 
is measured in units of . These potentials are related 
to each other: 

VaaiR,aa) = Vra{^R, aa)[l + 0{aaf] (25) 

= ^Var{R, aa) + AV;cif(a.)] [1 + Oiaa)% (26) 

where AV^eU = Vaa,se\f - £.Var,sc\f- While Vaa and Vra 
are equal at a given physical distance (up to lattice arte- 
facts) , in the case of Var a shift by an additive constant 
is expected since the self energies differ: 



AV;eif(a^) = 0.08214, 



(27) 



The numerical value has been obtained in lowest order 
perturbation theory for ^ = Co = 4. 

Following Ref. |£^ we use Eq . (25| ) to determine the 
renormalised anisotropic^ Eq. (|2£ ) can then be used as 



^ We ignore the possibility of mixing r and a coordinates 
within the "spatial" separation. 



Note that unlike Ref. 



our analysis is based on asymp- 



4 



an independent consistency check (modulo lattice arte- 
facts). In order to guarantee a consistent definition of ^ 
we can either consider the limit i? ^ oo or demand 
the matching to be performed at the same distance in 
terms of a measured correlation length. We follow the 
latter strategy and impose, 



(28) 



at i?^ K, R„i — (2/3)ra/aa where ro ~ 0.5 fm is the 
Sommer scale of Eq. ( [l8| ) . We restrict ourselves to on-axis 
separations and take = 2, 3, 4 for (3 — 5.8, 6.0 and 6.2, 
respectively. This choice is justified by our subsequent 
analysis where we find Rm = 2.05(2) and Rm = 2.03(2) 
on the 8^ and 12^ (3 = 5.8 lattices and Rm = 2.96(4) and 
Rjn = 4.07(4) at /3 = 6.0 and /3 = 6.2, respectively. The 
renormalised anisotropy ^ is then obtained from Eq. ( |2^ ) 
by interpolating the potential Vra according to three pa- 
rameter fits. 



R 



(29) 



The errors are obtained via the bootstrap procedure. 
The resulting anisotropics and fit ranges employed, R e 
[Rmin, Lt/2\, as wcll as the fit parameters in units of 
are displayed in Table [II. Tmin = ^min/flcr denotes the 
"temporal" separation from which onwards effective mass 
data [Eqs. ( |2^ ) and (p3|)] saturated into plateaus. 

The renormalised anisotropics ^ are also included in 
the last column off] Table fv|. In the third column of this 

are displayed while in 



table the one loop results [g6|6£^ 

the second last column mean field estimates |^,^ are 
shown, 



|3^r = l3V{Uar){Uaa). 



(30) 



The temporal and spatial average plaquettes {{Uar) and 
{Uaa)) are also included in the table. While the renor- 
malised anisotropics are underestimated by the one loop 
results by about 10 % they are overestimated by the mean 
field values by almost the same amount. Finally, in Ta- 
ble we compile the bare anisotropics ^o,4 = 4^o/C at 
which we should have simulated in order to achieve ^ = 4. 
In doing so we assume that our statistical uncertainties 
on ^ of order 1 % will dominate over variations of the 
ratios ^o/C under a change of by less than 2 %. 

After determining the anisotropics, the potential V^t 
is fitted to the parametrisation Eq. (29) for r > rL_^ 
r,„ = 2ro/3. The results are compiled in Table |V]| . 



The fit parameters e and K agree within errors with 
those determined from the data on Vra of Table HI 
while the Vq values tend to be somewhat smaller, in 
agreement with the expectation of Eq. (^). Note that 
the parametrisation is thought to be effective only and 
that the fit ranges employed for the two potentials dif- 
fer from each other. From the fit parameters, values 
ro/o,cr = \/(1.65^^i ^S}/^ ''^^ ^® extracted. These 
are displayed in Table p, along with the linear spatial 
lattice extent. Compared to the isotropic case, = 1; 
where ro/a = 3.64(5), 5.33(3) and 7.29(4) at 

P — 5.8,6.0 and 6.2, respectively, is somewhat in- 
creased while the temporal lattice spacing Ur is reduced 
The ratios rga^^^'^a 
exhibit that at ^ = 
tice spacings are about 15 % smaller than their isotropic 
counterparts, obtained at the same (3 values. 



= 4.36(4), 6.31(10) and 8.61(9) 
4 the geometrically averaged lat- 



TABLE III. Fits to V-rcr and determination of the renor- 
malised anisotropy ^. The fit parameters are displayed in 
units of a„. 








tmin 






Vo,T<J 




A"r<.C 


5.8 


8 


3 


6 


4.052(32) 


0.826(43) 


0.352(43) 


0.143(11) 


5.8 


12 


3 


6 


4.102(33) 


0.834(45) 


0.366(44) 


0.144(12) 


6.0 


12 


4 


5 


4.084(74) 


0.813(20) 


0.332(18) 


0.0697(49) 


6.2 


16 


4 


8 


3.957(51) 


0.792(18) 


0.356(28) 


0.0384(28) 



TABLE IV. Bare, one loop, and mean field es- 
timated ("tadpole improved") anisotropies versus the 
non-perturbatively determined renormalised anisotropy. 



/3 


Co 












5.8 


3.10 


3.470 


0.41557(3) 


0.807175(5) 


4.320 


4.076(25) 


6.0 


3.20 


3.573 


0.44540(2) 


0.820463(4) 


4.343 


4.084(74) 


6.2 


3.25 


3.619 


0.46924(1) 


0.829761(2) 


4.322 


3.957(51) 



TABLE V. Estimates of the bare anisotropies that corre- 
spond to ^ = 4. 



/3 


5.8 


6.0 


6.2 


Co, 4 


3.04(3) 


3.13(6) 


3.29(4) 



totic r ^ oo results rather than on pre-asymptotic finite T 
approximants to the potential. 

^Note that in this table we have averaged the results ob- 
tained on the 8^ and 12^ lattices at /3 = 5.8 that agree with 
each other within errors. 
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FIG. 2. The three potentials Wo-, Vo-o- and Vo-t in units of 
at /3 = 6.2, = 3.25. 
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FIG. 3. Differences between Vacr and ^Vo-t- 



TABLE VI. Fits to Vo-t. The fit parameters are displayed 
in units of tto-. 



/3 




^ mill / CIt 










5.8 


8 


13 


V5 


0.756(50) 


0.348(74) 


0.138(8) 


5.8 


12 


13 


V5 


0.732(55) 


0.298(80) 


0.145(9) 


6.0 


12 


10-11 


3 


0.745(39) 


0.354(80) 


0.0659(41) 


6.2 


16 


16-19 


4 


0.658(25) 


0.249(67) 


0.0376(17) 



Assigning the phenomenological value | |65| , |7l| , [7^ , [70t 
0.5 fm to To we find L^a^ « 1.3 fm on the small lattices 
and Lcacr w 2 fm on the 12'^ lattice at /? = 5.8. This 
means that ^/3L„acr/2 > 1.1 fm in all our simulations; 
along the R || (1,1,1) direction we are safe from the ef- 
fect of mirror charges up to distances bigger than one fm. 
Beyond this distance only representations with non-zero 
triality are protected by the centre symmetry ]70|] from 
direct finite size effects. 

In Fig. ^we display all three potentials in units of Oo- at 
(3 — 6.2. Note that the anisotropy has been determined 
by matching Vra- to Vo-o- at r = ia^- In addition to the 
data points two curves are included that correspond to 
the parameter values of Table III and Table VI from fits 
according to Eq. ( ^ to Vra for r > 2a^ and to Var for 
r > Aa^, respectively. The matched potentials Vaa and 
Vra follow the same curve. 

Up to lattice artefacts and the self energy shift AVscU, 
£,Var and Vaa , that both live along coarse grained lattice 
directions, should also agree with each other [Eq. (p6|)]. 
Indeed, as is demonstrated in Fig. |3[ the differences 
are compatible with constants of the order suggested 
by tree level perturbation theory for C = = 4, 
Eq. (|2^). Averaging the R > data points results in 
the values, 0.091(9), 0.096(18) and 0.124(11) for AFsoif 
at (3 = 5.8, 6.0 and 6.2, respectively (solid lines with er- 
ror bands). On the large lattice at /3 = 5.8 we obtain the 
value 0.081(9)a^^, in agreement with that above, from 
the smaller volume. These shifts of the self energies re- 
sult in reduced relative errors of {War) [Eq. (^] (and in 
increased errors of {Waa)), relative to the isotropic case. 



IV. THE POTENTIALS 

The potentials in non-fundamental representations are 
extracted in the same way as discussed above from fits 
to the corresponding smeared Wilson loops. These are 
obtained from the fundamental ones by use of Eqs. (0) 
- (|l7|). In the case of the fundamental Wilson loops, dis- 
cussed in Sec. [II, temporal links have been replaced by 



their thermal averages in the vicinity of the surrounding 
staples in order to reduce statistical fluctuations [|7llj73| 
(link integration |7^). Note that our use of Eqs. (|lO|) 
- ( p^ implies that we cannot thermally average funda- 
mental links in the construction of higher representation 
Wilson loops anymore. 

We determine the potentials from correlated exponen- 
tial fits to {WD,aT) data according to Eq. (||). The fit 
range in T is selected separately for each distance R and 
representation D, such that x^/^df < 1-5. In addi- 
tion, we demand the saturation of "effective masses", 
VD{T)aa = ^ln[{WD{T)y{WD{T + A))], into plateaus 



ranges T € [T, 



T 

mini max 



]t: 



for T > Tmin. In Table VII we display the resulting fit 



;hat have been selected by means 
of this procedure for the example of the point Vd (R) with 
R M ro/aa- In general, the interplay between statistical 



6 



errors and ground state overlaps was such that rmin only 
slightly varied with R. In the case of the fundamental 
potential we find values 2ro < tmin < 4.5 tq, depending 
on R and the parameter values we simulate at, while for 
D = 15s we find, 0.7 rg < imin < '"o- The corresponding 
esti mate d ground state overlaps cd are displayed in Ta- 
ble VII] . The overlaps decrease with increasing Casimir 
constant, lattice spacing or distance r. At r « rg the 
overlaps range from 0.62(3) in the worst case to 0.97(1) 
in the best case which quantifies the efficiency of our 
smearing algorithm. 

In Fig. ||we display the resulting potentials at /3 = 6.2. 
The curves correspond to the three parameter fit Eq. ( |29| ) 
to the fundamental potential, multiplied by the respec- 
tive ratio of Casimir factors, do of Table |. It is clear 
that the Casimir scaling hypothesis Eq. (^) works quite 
well on our /3 = 6.2 data for the investigated distances. 
In Figs, [s] - we display the ratios Vd {R, aa)/VF [R, a^) 
for our three lattice spacings. We did not attempt to sub- 
tract the self energy contributions [cf. Eqs. (|^) - (||)] in 
this comparison. At /3 = 5.8 and /3 = 6.0 we find the data 
to lie significantly below the corresponding Casimir ra- 
tios (horizontal lines). However, the deviations decrease 
rapidly as the lattice spacing is reduced. 



TABLE VII. Fit ranges tmin/ar - fmax/or at r ' 



■ ro. 



D 


P = 5.8 


P = 6.0 


P = 6.2 


3 


13 " 24 


10 - 36 


16 - 40 


8 


7-12 


9-18 


10 - 20 


6 


7-12 


8-18 


10 - 20 


15a 


6-11 


6-14 


8-15 


10 


5-10 


5-11 


6-12 


27 


4-8 


4-8 


5-9 


24 


4-8 


4-8 


5-9 


15s 


3-7 


3-7 


4-7 



TABLE VIII. Ground state overlaps c_d at r ~ ro. 



D 


P = 5.8 


P = 6.0 


P = 6.2 


3 


0.89 (1) 


0.96 (1) 


0.97 (1) 


8 


0.81 (1) 


0.90 (3) 


0.91 (2) 


6 


0.80 (2) 


0.90 (3) 


0.90 (3) 


15a 


0.68 (4) 


0.84 (3) 


0.88 (5) 


10 


0.62 (3) 


0.82 (2) 


0.89 (3) 


27 


0.67 (3) 


0.81 (2) 


0.88 (3) 


24 


0.65 (3) 


0.77 (2) 


0.90 (4) 


15s 


0.66 (3) 


0.79 (2) 


0.86 (3) 



CO 
CO 



> 




FIG. 4. The potentials for all measured representations, 
obtained a,t P — 6.2. Note that we did not subtract any self 
energy pieces but just rescaled the raw lattice data in units 
of ro. 
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FIG. 5. The potentials normalised to the fundamental po- 
tential Sit P — 5.8, in comparison to the expectations from 
Casimir scaling (horizontal lines). 
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FIG. 6. The same as Fig. at /? = 6.0. 
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FIG. 7. The same as Fig. | at /? = 6.2. 

Prior to a continuum limit extrapolation of the ratios 
we shall investigate finite size effects by comparing re- 
sults obtained on the 1.3 fm lattice with results from the 
2 fm lattice at /3 = 5.8. For the fundamental potential 
we already know from previous studies that for spatial 
extents, L^aa > 2ro, such effects are well below the 2 % 
level The situation is less clear for potentials in 

higher representations. In principle, the flux tube be- 
tween the sources could widen when the energy per unit 
length is increased and, therefore, higher representation 
potentials might be more susceptible to finite size effects. 

In Fig. H we compare the fundamental, octet and sextet 
potentials obtained on the 12'^ lattices at /? = 5.8 (full 
symbols) with those obtained on the 8'^ lattices (open 
symbols). Up to distances well beyond 2ro ~ Ga^- no 
statistically significant deviations are seen. In Fig. pi we 

^^'=i2^"^'-^-=«(r)-l, 



show the relative deviations, V, 



D 



between the potentials determined on the larger lattices 



and those measured on the smaller lattices for all the 
representations that we have investigated. Again, no 
systematic or statistically significant differences are de- 
tected. Up to r = 4 Oct ~ 1.4 ro this holds true on the 
1 % level for the fundamental potential and on the 3-5 % 
level for higher representation potentials. Beyond this 
distance the statistical errors start to explode. The same 
comparison has been performed for ratios of potentials. 
In this case the relative errors are slightly reduced due to 
correlation effects. However, no statistically significant 
tendencies were observed either. The relative statistical 
errors on the two lattices are of about the same size and 
comparable to those of the f3 = 6.0 and (3 = 6.2 simula- 
tions. Thus, we do not expect finite size effects to exceed 
the statistical errors on any of the simulated lattice vol- 
umes. 




FIG. 8. Comparison between the potentials in the lowest 
dimensional three representations between the La = 8 lattices 
(open symbols) and La — 12 lattices (full symbols) at /3 = 5.8 
in lattice units. 

We now attempt a continuum limit extrapolation of 
our data. We remark that in the limit ^ [Eqs. (|l|) 
- (j^)] the Casimir scaling of the diverging self ener- 
gies Vo.scuio-cr) automatically implies Casimir scaling of 
yD(R, Oct). However, with Vo.scif being a purely ultra 
violet quantity, this sort of Casimir scaling has little to 
do with non-perturbative aspects of the theory. As can 
be seen from Fig. ^, where the potentials vary by more 
than a factor two with the distance, even at our finest 
lattice resolution Vo.seif have not yet become the domi- 
nant contributions to the lattice potentials. In order to 
avoid the trivial Casimir scaling described above we will 
only study ratios of physical interaction energies. 



Roir) 



Voir) VD(R,a, 



Vr 



,sclf(Qa-) 



,self(a<T 



+ Oial)], 
(31) 



where r = Ra„ 



8 



0.05 


-0.05 

0.05 


-0.05 
0.1 


-0.1 
0.1 



-0.1 
0.1 



-0.1 
0.2 



-0.2 
0.2 



-0.2 
0.2 



-0.2 



- *H ,4 4 



, • • • »l i 



15a 



10 



15s 



- '27 

- • • • • • • i 















i 




- 24 

- • • • » • • • » « 




...».J 


* - 











1 



2 3 4 5 6 7 

r/a„ 

FIG. 9. Relative deviations yj^"=^2(r)/V^<'=**(r) - 1 of the 
potentials obtained on different volumes as a function of r/aa- 

We estimate the self energies in leading order perturba- 
tion theory, Eq. (^. Lattice perturbation theory is noto- 
rious for its bad convergence behaviour |^,^. However, 
in our anisotropic case the effective expansion parameter 
g'^£,Q^ is much smaller than in standard applications of 
lattice perturbation theory. In Fig. ^we have indeed seen 
that leading order perturbation theory predictions on the 
difference of two self energies agree within 30 % with nu- 
merical data. We estimate the self energies in two dif- 
ferent ways: (a) we use the bare lattice coupling and the 
renormalised anisotropy ^, (b) we mean field ("tadpole") 
improve ||6^,B^ both, coupling and anisotropy [Eq. (pO|)], 
^0ir9ir ~ ^0 y'^iUar)^^- The rcsults from the two meth- 



ods, shown in Table IX, differ by up to 60 % from each 
other. We will use the estimate (a) in our analysis but 
take the difference between (a) and (b) into account as 
a systematic error. While data at large distances are 
marginally affected by this uncertainty, the error bars at 
small distances are significantly increased. 



cc 




i i if — * Mi ^ 



15a 
6 



8 



1 1 : T i i. ■ -r 

6 

_ 5 
4 
3 
2 
1 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

r/ro 

FIG. 10. Continuum limit extrapolated ratios, Eq. (^l|), 
compared to the Casimir scaling hypothesis (horizontal lines) . 

After subtracting the (scaling violating) self energy 
contributions we determine the continuum extrapolated 
ratios Roir) by means of quadratic fits, Eq. (pi|). We 
perform these extrapolations for all the distances r that 
have been realised on our coarsest lattice (/3 = 5.8) in 
units of rg. On the finer lattices, we linearly interpo- 
late between the two lattice points that are closest to 
each given distance r, prior to the quadratic continuum 
limit fit. We find the data to be compatible with the 



g- 


10 


KI 





The numerical values are displayed in Tables ^ 
No statistically significant violations of Casimir 
scaling are found. Our accuracy is somewhat limited at 
short distances, due to the perturbative estimation of 
the self energies. The slope of the extrapolation in 
increases with the distance r as well as with the Casimir 
charge Cd', large masses VoiR^a^) > are more af- 
fected by lattice artefacts than small masses. This ob- 
servation also explains why the deviations from Casimir 
scaling at /3 = 5.8 and /3 = 6.0 increase with the distance 
(Figs. I - 1). 

TABLE IX. Perturbative estimates of the self energy 
Vp self (in units of a^). 



/3 


bare g^, renormalised ^ 


mean field gfr^^}^ 


5.8 


0.0855(05) 


0.13931 


6.0 


0.0825(15) 


0.12835 


6.2 


0.0824(11) 


0.12093 
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TABLE X. The continuum ratios, Roir) = VD{r)/VF{r). 
"exp." denotes the Casimir scaling expectation. 



TABLE XL The same as Table ^ for higher representa- 
tions. 



r/ro 


U — o 


U — D 


n ^ ^/-i 

U — lOd 


A QQ 
U.OO 


z.z4{^Uy ) 


O /I O 1 1 \ 

Z.4o(^iij 


o.y ( [zy ) 


n Ad 
U.4d 


z.z/ (Uoj 


o c: o / A /I \ 
z.0o(U4j 


o OA / 1 n\ 

o.yy(iuj 


U.OD 


Z.21(U4j 


o c: a/a /I \ 
z.0U{^U4j 


o / 1 c:\ 
O.OD(l0 ) 




O OO / A /I \ 

Z.2o(U4j 


z.4o{^Uoj 


O A^ / 1 1 \ 

o.yo^ii ) 


A 'VO 


O O /I / AO^ 

z.24(U2j 


z.oU^Uoj 


o at/AO \ 

O.y (^(^Uoj 


A OA 




O /I K / A /I \ 

2.4o(^U4j 


o A c: / 1 1 \ 

o.yo(^ii j 


A AO 


z.z4((Joj 


O /I /AC \ 

z.4o(Uoj 


A AA/A0\ 

4.UU(Uoj 


U.yo 


Z.Zo(^U4j 


KA/np;\ 
Z.OKjyJO ) 


Q QV/ 1 0\ 


A AO 


z.zz(^U4 j 


z.4y(U4j 


Q Q / 1 1 ^ 
O.O0(^ii ) 


1 1 Q 

i.io 


O 0/1 /'A/1 \ 

z.z4(^U4 j 


z.4o(U0 ) 


Q QQ / 1 0^ 

o.oo(^iz ) 


1 QA 

i.oU 


Z.Z / (^Uo ) 


z.4y(^u / ) 


4.11(^iD J 


1 QQ 

l.of5 


Z.Zo(^U4j 


/1Q/'^7^ 
Z.4o(^U / ) 


A AC / 1 


1 AC 

i.4o 


O O /I / a /I \ 

Z.Z4(U4j 


o c: c; /AC \ 

z.ool^Uoj 


O O^i ( 1 o\ 
O.OD(loJ 


1 c:a 

i.oy 


z.z4(^U0 j 


O /I A/ AC\ 

z.4y(^Uo j 


A 11/1 A^ 

4.ii(iy ) 


i.Oo 


z.zy^uo ) 


O fiA/ AC\ 

z.oUl^Uo j 


A AA/ 1 Q^ 

4.uy(io ) 


i.oy 


Z.Zo(^Ud ) 


Z.DlyW ) 


Q Q/1 / 1 




Z.ool^UD j 


9 p:c/aq\ 
z.ooi^uy ) 




1.95 


2.26(10) 


2.67(12) 


4.10(19) 


1.95 


2.27(09) 


2.66(11) 


4.10(18) 


2.18 


2.30(11) 


2.11(16) 


3.67(26) 


2.25 


1.97(16) 


2.03(20) 


3.41(43) 


2.30 


2.26(14) 


2.45(17) 


3.72(31) 


2.39 


2.29(14) 


2.46(16) 




exp. 


2.25 


2.5 


4 





_D = 10 


D = 27 


D = 24 


D = 15s 


0.33 


4.45(37) 


6.23(65) 


5.98(57) 


6.84(83) 


0.46 


4 39Q8~) 


6 10f33~) 


5 94('22'l 


6 95('27'l 


0.56 


4 38Q6~) 


5 96f2Q~) 


5 82('23'l 


6 87('23'l 


0.65 


4.37(13) 


6.10(16) 


5.86(14) 


6.86(19) 


0.73 


4.45(11) 


6.21(15) 


5.97(14) 


6.90(17) 


0.80 


4.43(12) 


6 18Q6~) 


5 92(14") 


6 89(17") 


0.92 


4.43(10) 


6.17(16) 


5.95(14) 


6.81(21) 


0.98 


4.34(14) 


6.28(14) 


5.95(15) 


7.00(13) 


0.98 


4.37(14) 


6.21(12) 


5.91(14) 


6.95(12) 


1.13 


4.48(10) 


6.06(17) 


5.82(19) 


6.71(19) 


1.30 


4.46(11) 


5 76f24") 


5 52(241 


6 99(16") 


1.38 


4 5nnn") 


fi 00(^21 ") 


5 99f23"l 


7.08(14) 


1.45 


4.44(10) 


6.21(16) 


5.96(18) 


7.19(12) 


1.59 


4.55(12) 


6.50(23) 


6.52(23) 


7.27(15) 


1.63 


4.48(13) 


6.02(21) 


6.05(29) 


7.23(14) 


1.69 


4.50(16) 


6.39(31) 


6.19(28) 


7.33(17) 


1.84 


4.79(14) 


6.53(27) 


6.09(22) 


7.51(16) 


1.95 


4.74(20) 


6.51(33) 


6.21(30) 


7.16(48) 


1.95 


4.71(19) 


6.54(32) 


6.16(29) 


7.33(35) 


2.18 


4.40(36) 


5.93(53) 


5.89(43) 


7.19(57) 


2.25 


3.12(50) 


5.67(71) 


5.99(58) 




2.30 


3.82(44) 


6.37(76) 


5.97(53) 




2.39 


3.62(53) 








exp. 


4.5 


6 


6.25 


7 
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V. DISCUSSION 

We have confirmed that violations of the Casimir scal- 
ing hypothesis, 



(32) 



are below the 5 % level for distances r < 2ro « 1 fm 
in the continuum limit of four dimensional SU (3) gauge 
theory for all representations with Cd l£ 7- This find- 
ing rules out many models of non-perturbative QCD and 
imposes serious restrictions onto others. For instance in 
a bag model calculation scaling of string tensions with 
the square root of the respective Casimir ratio has been 
obtained [75|] and instanton liquid calculations result in 
ratios between potentials in different representations that 
are smaller than the Casimir ratios too . 

Another possibility would have been scaling propor- 
tional to the number of fundamental flux tubes em- 
bedded into the higher representation vortex [p -\- q in 
SU{3)] ^9 41 1, which happens to coincide with Casimir 
scaling in the large N limit of SU{N). This picture is 
supported by the finding that the SU{N) vacuum seems 
to act like a type I superconductor ||77| , |78| , i.e. flux tubes 
repel each other. However, this scenario is also excluded 
by the present study. Furthermore, serious restrictions 
onto most of the remaining models are imposed (see for 
instance Ref. |5^). It is particularly disappointing that 
neither centre vortex models ||47|-|49[| nor the dual super- 
conductor scenario or string models [0 seem 
to offer any explanation why the numerical data so closely 
resemble the Casimir ratios. Certainly, it is worthwhile 
to dedicate more theoretical effort to this fundamental 
phenomenon. 

We have not discussed "string breaking" so far. While 
the fundamental potential in pure gauge theories lin- 
early rises ad infinitum, the adjoint potential will be 
screened by gluons and, at sufficiently large distances, 
decay into two disjoint gluelumps This string 

breaking has indeed been confirmed in numerical stud- 
ies ||,^,|^,||,|||. Therefore, strictly speaking, the ad- 
joint string tension is zero. In fact, all charges in higher 
than the fundamental representation will be at least par- 
tially screened by the background gluons. For instance, 
6(8)8 = 24® 15a* ©6© 3*: in interacting with the glue, 
the sextet potential obtains a fundamental component. A 
simple rule, related to the centre of the group, is refiected 
in Eqs. (10) - (|l7|): wherever zP~'' = 1 (zero triality), 
the source will be reduced into a singlet component at 
large distances while, wherever z^"'' = z (or z*), it will 
be screened, up to a residual (anti-)triplet component, 
i.e. one can easily read off the asymptotic string tension 
(either zero or the fundamental string tension) from the 
third column of Table |, rather than having to multiply 
and reduce representations. As a result, the self-adjoint 
representations, 8 and 27, as well as the representation. 



10, will be completely screened while in all other repre- 
sentations with p -\- q < 4 a residual fundamental compo- 
nent survives. The same argument, applied to SU(2), re- 
sults in the prediction that all odd-dimensional (bosonic) 
representations arc completely screened while all even- 
dimensional (fermionic) representations will tend towards 
the fundamental string tension at large distances. 

One expects this sort of string breaking and fiattening 
of the potential to occur at distances larger than about 
2.4 rp |3^. Obviously, once the string is broken Casimir 
scaling is violated. It is certainly interesting to investi- 
gate what happens around the string breaking distance. 
However, this requires lattice volumes exceeding those 
used in the present study as well as additional operators 
that are designed for an optimal overlap with the respec- 
tive broken string states [p2[. 
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